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^ ! Abstract 

<n : 

I | Recent progress has been experienced in the field of hypernuclear weak de- 

cay, especially concerning the ratio of the neutron- to proton-induced A non- 

^ ' mesonic decay rates, T n /T p . Theoretical analyses of nucleon coincidence data 

■ have been performed in a finite nucleus framework. They led to the extraction 

of T n /T p values in agreement with pure theoretical estimates, thus providing 

\Q ', an evidence for the solution of a longstanding puzzle. Here we present an alter- 

native approach to the problem, based on a nuclear matter formalism extended 

i -^h ' to finite nuclei via the local density approximation. The work is motivated by 

1 { ■ the exigence to make the determination of T n /T p from data less model depen- 

dent. One-meson-exchange potentials are used for describing both the one- 
and two-nucleon induced decays, AN — > nN and ANN — > nNN. For the 
latter, treated within a microscopic approach, the channels Ann — > nnn and 

J^j ' App — > npp are included in addition to the mode Anp — > nnp already consid- 

5_i ■ ered, in a phenomenological way, in previous studies. The propagation of the 

final nucleons in the residual nucleus is simulated by an intranuclear cascade 
code. We evaluate single and double coincidence nucleon spectra for the non- 
mesonic decay of ^ 2 C. Through the comparison of our predictions with KEK 
coincidence data we determine T n /T p = 0.43 ± 0.10 for this hypernucleus, 
confirming previous finite nucleus analyses. The use of a high nucleon kinetic 
energy detection threshold such as 60 MeV makes the contribution of two- 
nucleon induced channels quite small, but final state interaction effects are 
still important when extracting r n /r p from measured single nucleon distribu- 
tions. Coincidence spectra suffer less from these effects even for a moderate 
detection threshold such as 30 MeV. In any case, final state interactions have 
to be considered for meaningful determinations of r n /r p . 
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I. INTRODUCTION 



Diversified efforts have been devoted to the study of hypernuclear weak decay in the 
latest years. Theoretical reviews on the subject can be found in Refs. [1,2] and recent, 
related experiments in Refs. [3-10]. Forthcoming data are expected from FINUDA [11], 
while new experiments are planned at J-PARC [12] and HypHI [13]. Strong evidences 
for a solution of the longstanding problem on the ratio r„/r p between the widths for the 
weak processes An — > nn and Ap — > np originated from theoretical analyses [14,15] of 
KEK nucleon coincidence data [7-10]. In the authors opinion, this puzzle was due to the 
non-trivial interpretation of experimental data, which required a careful analysis of nuclear 
medium effects on the weak decay nucleons, rather than to a poor understanding of the 
weak decay mechanism. Indeed, the lately extracted T n /T p values [14,15] turned out to 
be in agreement with the previous, pure theoretical estimates of Refs. [16-20] obtained by 
using one-meson-exchange potentials to describe the one-nucleon induced, AN — > nN weak 
transitions. 

Nevertheless, further theoretical and experimental work is desirable in order to confirm 
the previously mentioned evidence in favor of a solution of the r n /r p puzzle. Indeed, on the 
one hand such an evidence relies on particular theoretical descriptions of both the weak decay 
mechanism and the subsequent propagation of the produced nucleons within the residual 
nucleus. In this direction, the use of alternative weak decay or/and intranuclear cascade 
models is of interest. On the other hand, one has to consider that another problem of the 
field is still unsolved: it concerns the asymmetry of the protons emitted in the non-mesonic 
decay of polarized hypernuclei, measured to be not far from zero in recent experiments 
[21,22] while a large negative number is predicted by the theoretical models [17,23-26]. 
Recently, a strong effect of nucleon final state interactions (FSI) was pointed out [23] without, 
however, bringing new hints for a possible solution of the asymmetry puzzle. The connections 
existing among the weak decay observables [r n /r p , T^m = T n + T p + T2, T2 being the width 
for two-nucleon induced decays, ANN — > nNN, and the asymmetry parameters] and the 
question concerning the validity of the AI = 1/2 isospin rule in the non-mesonic decay 
is another important issue which deserves future investigations in the prospect of a better 
understanding of baryon-baryon weak interactions. 

On the experimental side, very recent measurements of single- [3-6] and double- 
coincidence [7-10] nucleon spectra from the non-mesonic hypernuclear decay were reported 
— with accuracies largely improved with respect to the ones at disposal in previous exper- 
iments [27-31] — in forms that suggest suitable comparisons with theory. Some of these 
experiments have somehow managed to derive values of r n /r p . They are discussed in the 
following paragraphs. 

The T n /T p ratio was obtained by KEK-E307 [3,4] for ^ 2 C, ^ 8 Si and AFe hypernuclei from 
single-proton kinetic energy spectra measurements and by making use of the intranuclear 
cascade code of Ref. [32] (based on the polarization propagator formalism of Ref. [33]) 
to simulate the spectra of the nucleons emitted by the considered hypernuclei. For \ 2 C, 
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ryr p = o.87 ± 0.23 (o.601 oil!) was obtained [4] by neglecting the two-nucleon induced 
decay mechanism (for T2/(r n + T p ) = 0.35, with T2 = T(Anp — > np)). As the very same 
authors of Refs. [3,4] noted in Ref. [5], these determinations of the ratio may be affected by 
the fact that in the experiment the neutron-induced decay width was estimated indirectly, 
from the proton measurement, using the relation r„ = I\ — T p — F n - — F n o, which neglects 
two-nucleon stimulated non-mesonic decays. This method also required the measurement 
of the total decay width T T as well as the decay rates for the mesonic channels A — > ir~p 
(T w -) and A — > ir°n (r^o, for which previous data from Ref. [28] were used in the analysis of 
Ref. [4]). Moreover, the severe energy losses suffered by protons inside the (thick) target and 
detector materials and the consequently high kinetic energy threshold (about 40 MeV) for 
proton detection in KEK-E307 did not permit an easy reconstruction of the proton spectrum 
emitted by the nucleus, which is essential for the indirect evaluation of T n . Theoretical input 
about nucleon rescattering in the residual nucleus [32] was indeed necessary to supply to 
this problem. As a consequence, in Ref. [5] the hypothesis was advanced that T n (T p ) might 
be overestimated (underestimated) in the analysis of Ref. [4] because of an underestimation 
in the number of emitted protons. 

A controversial determination of the ratio from KEK-E369 data, based on non- 
demonstrated, delicate hypotheses and theoretical input (again from Ref. [32]), was reported 
in Ref. [5]. In this experiment, direct measurements of single-neutron kinetic energy spectra 
were performed (with a 10 MeV threshold) for ^ 2 C and ^ 9 Y; once analyzed together with 
the single-proton spectra of Refs. [3,4], a ratio T n /T p = 0.51 ± 0.15 for ^ 2 C was derived by 
neglecting the two-nucleon induced decay channel. We shall comment on the reliability of 
the extraction method used for such a determination in Section III B. 

To overcome the difficulties of the discussed KEK experiments, both single-neutron and 
single-proton energy spectra were measured simultaneously by KEK-E462 for ^He and 
KEK-E508 for ^ 2 C [6]. From these measurements, the authors concluded that T n /T p ~ 
(N n /N p — l)/2 ~ 0.5 for both \Re and \ 2 C, N n (N p ) being the total number of neutrons 
(protons) with kinetic energies T/v above 60 MeV. However, one has to note that the previous 
approximate relation between T n /T p and N n /N p is only valid when FSI and two-nucleon 
induced decay effects can be neglected. The predictions of Ref. [15] and the results presented 
in Sections III A and III B prove that FSI are not negligible even when a high detection 
threshold such as T^ 1 = 60 MeV is used. 

In the experiments KEK-E462 (^He) and KEK-E508 (\ 2 C), nucleon-nucleon coincidence 
spectra were also measured [7-10]. Quite clean angular and energy correlations between 
neutron-neutron and neutron-proton emitted pairs (i.e., back-to-back kinematics and T/v x + 
Tjv 2 — 155 MeV) were observed, thus representing the first direct experimental evidence 
of the existence of the two-body decays An — * nn and Ap — > np. The ratio, N nn /N np , 
between the numbers of emitted neutron-neutron and neutron-proton pairs was measured 
to be around 0.5 for both |He and ^ 2 C after applying the angular and energy restrictions: 
cos9nn < —0.8 and T N > 30 MeV. The authors of Ref. [8,9] concluded that, under these 
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constraints, T n /T p ~ N nn /N np ~ 0.5 on the basis of a supposed cancellation of FSI and 
two-nucleon stimulated decays effects. In a very recent work [10], the result T n /T p = 0.51 ± 
0.13 ± 0.04 was deduced by the KEK collaboration for ^ 2 C after correcting for FSI effects 
by making use of the number of detected proton-proton pairs in addition to measurements 
of N nn and N np . A rather schematic method, not completely reliable in our opinion, quite 
in line with the one used in Ref. [5], was applied to determine T n /T p . We shall discuss the 
effect of FSI and two-nucleon induced decays in the extraction of r n /r p from measurements 
of N nn /N np in Section HIE. 

Despite this recent experimental progress, improved and/or independent measurements 
are awaited for a really complete understanding of the AiV — > nN reaction in nuclei. On this 
respect, the observation of the weak decay of neutron and proton rich hypernuclei [11,13] 
would also be source of new information. 

On the theoretical side, Refs. [14,15] presented an extensive study of single- and double- 
coincidence nucleon spectra for the non-mesonic decay of ^He and ^ 2 C hypernuclei. A one- 
meson-exchange (OME) model was used for the AA" — > nN transition in a finite nucleus 
framework. The two-nucleon induced decay channel Anp — > nnp was taken into account via 
the polarization propagator method in the local density approximation of Refs. [33,34]. The 
intranuclear cascade code of Ref. [32] was used to simulate the nucleon propagation inside the 
residual nucleus. Comparison with KEK-E462 and KEK-E508 coincidence data [7-9] lead 
to the determination of T n /T p values around 0.3-0.4 for both |He and ^ 2 C. The relationship 
between T n /T p and the observable ratio between neutron-neutron and neutron-proton pairs, 
N nn /N np , was established for the first time in Ref. [15]. It was shown that FSI and two- 
nucleon induced decays significantly affect the extraction of T n /T p from coincidence data 
even when favorable energy and angular correlation restrictions are imposed on the observed 
nucleon pairs. Coincidence measurements of course make the determination of r„/r p easier 
and cleaner, since FSI and two-nucleon induced decay effects are reduced with respect to 
the ones emerging in analyses of single nucleon observables such as N n /N p . Nevertheless, 
contrary to what claimed in Refs. [7-9], they do not permit an exclusive identification of the 
non-mesonic decay channels as neutron- and proton-induced. Thus, r„/r p ^ N nn /N np (see 
figures 11 and 12 of Ref. [15] and the present discussion in Section HIE), indicating that to 
determine T n /T p one has to rely on experimental as well as theoretical coincidence spectra. 

With the purpose of making the extraction of T n /T p less model dependent, in the present 
paper we make use of an alternative framework for hypernuclear decay: it consists of a nu- 
clear matter formalism [19,20] extended to finite nuclei via the local density approximation, 
with the same weak OME transition potential (containing ir, p, K, K*, u and n exchange) 
of Ref. [17], in addition to the Monte Carlo intranuclear cascade code of Ref. [32]. At vari- 
ance with previous analyses, in the present paper we follow a microscopic approach for the 
two-nucleon induced decay, also including the channels Ann — > nnn and App — > npp besides 
the standard mode Anp — > nnp. Most of the results shown are for the intermediate-mass 
hypernucleus ]^C, but we also present some results for a heavier hypernucleus, ^ 9 Y. Both 
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cases can be well described within a local density approximation. 

The paper is organized in the following way. The weak decay models employed in the 
calculation are outlined in Section II A. In Section II B we give a few details of the in- 
tranuclear cascade simulation. Numerical results for single and double coincidence nucleon 
distributions are presented and compared with data in Section III. The contribution of the 
two-nucleon induced decay channels is analyzed with special regard. Finally, in Section IV 
we draw our conclusions. 



II. MODELS 

A. Weak decay 

The weak decay transitions as well as the distributions of the weak decay nucleons are 
obtained by means of a many-body description of the A self-energy in nuclear matter, based 
on the polarization propagator method (PPM) originally proposed by Oset and Salcedo [35]. 
The local density approximation (LDA) is employed to extend the calculation to finite nuclei. 
This approach was previously established [19,20] to evaluate the non-mesonic decay widths of 
\ 2 C There is a major value of this model which by itself adds sufficient novelty to the present 
work with respect to previous analyses [14,15] of the nucleon spectra from hypernuclear non- 
mesonic decay. We must indeed emphasize that the present approach is more microscopic 
(i.e., less phenomenological) than the previous ones. Indeed, in addition to the generally 
considered Anp nnp channel, also the Ann — > nnn and App — > npp mechanisms are now 
evaluated. The present approach is then more in line with the functional approach — in the 
framework of the bosonic loop expansion — used in Ref. [36], with the difference that a clear 
separation into one- and two-nucleon induced channels is now possible. 



1. One-nucleon induced decay 

We start with a brief summary concerning the evaluation of the one-nucleon (IN) in- 
duced decay widths T n and T p . Details can be found in Refs. [19,35]. It is convenient to work 
first with the partial decay width TtAtN—^t N it N n 

(p A , k F ), where p A is the A energy-momentum, 
kp is the Fermi momentum of nuclear matter and the Vs represent the isospin projections 
of the baryons. Using the standard Goldstone rules for diagrams in nuclear matter, it is 
straightforward to write for the partial decay width: 

(27T) 1 J (27T^ Gpart ^ A ~ ^ ^tfe) G ho l c (p 2 ~ q) 

4E (|(7 P i7 p2 |^" Af (g)| 7 A7, 2 )| 2 (2.1) 
-(7 P i7 P 2|U A ^" Af (g)| 7 A7,2)*(7 P 27 P i|^ A ^" Ar (pA -p 2 )| 7 A7/ l2 )) , 
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FIG. 1. Direct and exchange A self-energy diagrams corresponding to the one-nucleon induced 
decay channel in nuclear matter. 



where Pi (hi) stands for the energy-momentum of particles (holes). The meaning of each 
and hi is shown in Fig. 1, where we have drawn the A self-energy diagrams from which the 
relevant transition amplitude is obtained. From the energy-momentum conservation in each 
vertex we have p\ = p\ — q and h 2 = P2 — Q- For simplicity, 7, represents the spin and isospin 
projections of particle i. The one-meson-exchange weak transition potential, v AN ~* nN , 
takes into account the complete pseudoscalar and vector meson octets (tt, 77, K, p, u>, K*) 
through the parameterization of Ref. [17]. The summation in Eq. (2.1) runs over all spins 
and isospins of the weak transition potential. Note that the first term in the r.h.s. of this 
equation is the usually called direct contribution, while the second one is the exchange term. 
The particle and hole propagators are, respectively: 



\p\ ~ fcf 
p - E N (p) -V N + ie 
6{k F -\h\) 



G P M = _ f; ■ - , (2-2) 

Ghole(h) 



ho — E N (h) — Vn — is ' 

where E N (p) = m N +p 2 /2m^ is the nucleon total free energy and VJv the nucleon binding 
energy. 

By integrating Eq. (2.1) over p A one obtains the /c^-dependent width: 

r *At M -fpifp2 ( k F) = J dp A T tAth2 ^ tpltp2 (p A , k F ) |^a(Pa) I 2 , (2.3) 

where for the A wave-function, if) A (p A ), we take the lsi/2 wave-function of a harmonic 
oscillator and the A energy, (pa)o, is taken as (pa)o = n^A +pi/2mA + V\, where Va is the 
experimental A binding energy. To evaluate the decay width for a particular nucleus, one 
uses either an effective Fermi momentum or the LDA [35]. In the last case, kp is spatially 
dependent and the transition rate reads 

(Mr))|^ A (r)| 2 , (2-4) 
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where ipA(r) is the Fourier transform of ip\(p A ). 

Finally, by summing up over the isospin of baryons one obtains: 




(2.5) 



np 



and the lA-induced decay width is: 



r\ = r 



n 




(2.6) 



We have to stress that exchange terms are very important for an accurate evaluation of 
Y\ [19]. The present work calculates exactly the exchange terms in nuclear matter. The 
Random Phase Approximation (RPA) corrections have not been taken into account, because 
a precise evaluation of RPA-exchange terms is a rather involved task, which goes beyond 
the scope of the present contribution. 

At this point we have to mention that the models of Refs. [17,19,20] have been re- 
cently improved. On the one hand, a more realistic A wave-function, obtained in terms of 
the experimental hyperon binding energy, has been considered in the LDA calculation of 
Refs. [19,20], leading to a reduction of about 25% on the lA-induced decay rates, T n and 
T p , and of 30-40% on the 2AMnduced ones, T nn , T np and T pp , with respect to the original 
results. On the other hand, a numerically more accurate evaluation of the distorted final 
state wave-functions has been implemented in the finite nucleus approach of Ref. [17] and 
the rates increase by about 15%. In both cases the ratio T n /T p is essentially left unchanged. 
The updated predictions, used in the present work, are given in Table I. Fortunately, the 
spectra discussed in Refs. [14,15], being normalized per non-mesonic weak decay, are not 
affected by the corrections of Ref. [17]. Note that the non-mesonic rate T^m of the finite 
nucleus calculation is underestimated by about 25% because it lacks the contribution of the 
2A-stimulated mechanism. On the other hand, the LDA result is overestimated because 
two-body final state interactions between the emitted nucleons were not considered. As 
shown in Ref. [17], the omission of this effect leads to rates larger by as much as a factor 
of two. However, this factor disappears in the normalized spectra discussed in the present 
work and used for the important purpose of extracting r n /r p from measured distributions. 

For a comparison with experiment, in Table I we also report data for T NM and T n /T p . 
Only experimental determinations of T n /T p obtained from single-nucleon measurements are 
quoted. Large experimental errors affect all but the most recent data, especially for r n /r p . 
Almost all Y n /Y p data appear to strongly overestimate any theoretical prediction found in 
the literature. 



In the quasi-deuteron approximation, the 2AMnduced decay mode turns out to be dom- 
inated by the process Anp — > nnp: the meson emitted in the A vertex is mainly absorbed 



2. Two-nucleon induced decay 
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TABLE I. Non-mesonic weak decay rates (in units of the free A decay width) predicted for 
\ 2 C by the updated finite nucleus approach (OMEa and OMEf calculations) of Ref. [17] and LDA 
model of Refs. [19,20]. 



Ref. r n 


r„ 

V 


r r r 

- 1 - nn np 1 pp 




1 KTA/f 
IN 1V1 


Finite Nucleus, OMEa 0.190 


0.625 




0.303 


0.815 


Finite Nucleus, OMEf 0.173 


0.484 




0.356 


0.657 


LDA 0.267 


0.936 


0.017 0.238 0.062 


0.285 


1.521 


KEK-E508 [8] 








0.953 ± 0.032 


KEK-E369 [5] 






0.51 ±0.15 




KEK-E307 [4] 






0.87 ±0.09 ±0.21 


0.828 ± 0.056 ± 0.066 


KEK [30] 






1.87 ± 0.59tioo 


0.89 ±0.15 ±0.03 


BNL [29] 






1 qq + 1.12 

1 -°°-0.81 


1.14 ±0.20 



by an isoscalar neutron-proton correlated pair. However, the isospin quantum number also 
allows two others mechanisms, Ann — > nnn and App — > npp, which then contribute to the 
decay rate r 2 . Here we briefly outline our main expressions needed for the evaluation of this 
rate, which include all the just mentioned processes within a nuclear matter framework. A 
complete description of this decay channel can be found in Ref. [20] . 

At variance with r\, we limit ourselves to direct terms. This is because T 2 is originated 
from ground state correlations (GSC). It was shown in Ref. [37] that, in the case of electron 
scattering, the exchange terms for the two particle-two hole polarization propagator can be 
neglected. In fact, in that work it was concluded that exchange terms are not relevant in 
the graphs originated from GSC. 

Furthermore, it is a good approximation [38] to consider only contributions where the 
two weak transition potentials y hN ^ nN (which again include the exchange of ir, p, K, K*, 
rj and iq mesons [17]) of the A self-energy are attached to the same bubble. There are thus 
three different self-energy contributions, which we denote by pp, ph and hh. In the first one, 
the two V kN ^ nN are attached to the same particle (see Fig. 2). In the ph contribution one 
ykN^nN j g connec t e d to a particle and the other one to a hole. Finally, the two potentials 
are attached to the same hole for the hh part. To illustrate the method, we present the pp 
contribution to T 2 . From Goldstone rules, the pp-partial decay width is: 

f pp i u \ _o^ f d4p2 f d4p3 f d4q f d4q ' (0 7\ 

Lt A t h2 t h3 ^t pl t p2 t p3 {PA,kF) - ^mj (2 ^ )4 J (2 ^ )4 J ( ^ )4 J (2 ^ 4 

G Part(PA - q) G part {p 2 ) G 2 paTt {p 2 - q) G part (p 3 ) 

xGhoicfe -q + q') Choicfe - q') 



X 



(7pi7 p2 |V A ^" JV (g)| 7 A7 P '2) " (l P >2l P z\V NN U)\lh2l h z 



rNN i 



The meaning of each energy-momentum pi and hi is shown in Fig. 2. From energy- 
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FIG. 2. pp-part of the two particle-two hole contribution to the A self-energy in nuclear matter. 



momentum conservation in each vertex we have: p\ = p\ — q, p' 2 = P2 — q, ^2 = P2 — q + q' 
and h 3 = p 3 — q' . Again, r ) i stands for the spin and isospin projections of particle i. The 
nuclear residual interaction V NN is modeled by a Bonn potential [39] in the framework of 
the parameterization presented in Ref. [40], which contains the exchange of ir, p, a and uj 
mesons and neglects the small n and 5-mesons contributions. The summation in Eq. (2.7) 
runs over all spins and over t p / 2 . 

After performing the LDA as described in the previous subsection, the summation over 
the isospin of baryons leads to: 

fpp = fpp (9 %) 

nn Ann^nnn 5 V / 

ppp _ -pPP _|_ -pPP | -pPP _|_ -pPP | -pPP _j_ pPP 

np Anp— *nnp ' Apn^nnp ' Anp^npn ' Apn— *npn ' Anp^pnn ' Apn^pnn ? 

ppp _ -pPP | -pPP 

pp App^npp ' App^pnp • 

Thus, once the ph and hh channels are included, the 2AMnduced decay width is obtained 
as: 



r2 — r«ri ~\~ r?ip ~\~ rpp > (2-9) 

where (ij = nn, np or pp): 

r„ r;r • \t ■ vf. (2-io) 

Explicit expressions for all contributions can be found in Ref. [20] . 

The results obtained for the 2A r -stimulated decay widths are given in Table I. As ex- 
pected, the dominant contribution to T 2 originates from the np-induced decay, Anp — * nnp. 
Nevertheless, the other 2iV-induced channels cannot be neglected. The prediction T 2 /ri = 
0.26 is in agreement with the phenomenological estimate I^/Ti = 0.25 of Refs. [14,15]. 
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B. Intranuclear cascade 



After the A interacts with one (two) nucleons, the two (three) produced nucleons will 
interact with other nucleons in their way out of the nucleus. This process, which will 
generate secondary nucleons, is accounted for by the intranuclear cascade model described in 
Ref. [32]. This model considers a semiclassical propagation of primary (i.e., weak decay) and 
secondary nucleons. Nucleons move along classical, straight trajectories between collision 
points and under the influence of a local (i.e., independent ) mean potential, Vn(R) = 
—kp N (R) 2 /2rriN, where kp N (R) = [37r 2 p7v(-R)] 1//3 is the local nucleon (N = n, p) Fermi 
momentum corresponding to the nucleon density pn(R). Propagating nucleons also collide 
with the nucleons of the medium according to free space nucleon-nucleon cross sections 
properly corrected to take into account the Pauli blocking and Fermi motion effects. For 
more details of the code we refer to Ref. [32]. 

The description in terms of semiclassical nucleon propagation is justified when the nu- 
cleon wavelength X is much smaller than the average distance between nucleons, d, and the 
range of the nucleon-nucleon interaction, r : A << r < d [41]. For a 30 MeV kinetic energy 
nucleon, A = 0.8 fm, which has to be compared with d ~ 2 fm in \ 2 C and r ~ 1.4 fm. The 
spectra for nucleon kinetic energies < 30 MeV could thus show unrealistic behaviors. For 
this reason, all the comparison with data will intervene by introducing a kinetic energy cut 
of at least 30 MeV. 

III. RESULTS 
A. Single nucleon spectra 

In Figure 3 a comparison is shown between the present LDA calculation and the previous 
finite nucleus evaluation of Refs. [14,15]. The number of primary protons emitted in the IN- 
induced non-mesonic weak decay of \ 2 C is given as a function of the proton kinetic energy. 
In order to make the comparison model independent, the spectra are normalized per IN- 
induced decay assuming T n /T p = (r„/r p ) LDA = 0.285, which is equivalent to normalize 
per the corresponding proton-induced decay rate. We observe that the Full Width a Half 
Maximum (FWHM) of the LDA distribution is 15-20 MeV larger than the one of the finite 
nucleus spectrum. This is essentially due to the more pronounced Fermi motion effects in 
the nuclear matter calculation. Consistently with this expectation, in the finite nucleus 
calculation of Refs. [14,15] we found that the primary nucleon distributions for \ 2 C were 
slightly broader than the ones for |He. 

In Figure 4 we show the single-proton kinetic energy spectra for the non-mesonic weak 
decay of ^ 2 C once 2AMnduced decays and FSI effects are included. All spectra are normal- 
ized per non-mesonic weak decay. The continuous line refers to the present LDA result. The 
dashed line has been obtained with the LDA of Ref. [32] , where the weak transition potential 
was described by a correlated pion exchange. The old LDA spectrum corresponds to the 
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- - Finite Nucleus 

— LDA 




160 



T(MeV) 



FIG. 3. Kinetic energy spectra of primary protons from the lA-induced non-mesonic weak 
decay of jfC. The continuous line refers to the present LDA calculation, the dashed one to the finite 
nucleus evaluation of Refs. [14,15]. Both curves are normalized per lAMnduced decay assuming 

r„/r p = (r„/r p ) LDA = 0.285. 



same T n /T p of the present LDA, T n /T p = (r n /r p ) LDA = 0.285. Good agreement is obtained 
between the two LDA calculations despite the different weak decay models employed. The 
dot-dashed line has been taken from the finite nucleus evaluation of Refs. [14,15] fixing 
Tn/r p = (r«/rp) LDA = 0.285. Again, the differences between LDA and finite nucleus esti- 
mates are mainly due to the different phase spaces in the two cases. As it is apparent from 
Figure 4, all the theoretical spectra are in strong disagreement with KEK-E508 data [6]. 
Not even a calculation enforcing a large value of T n /T p such as 1 in the old LDA (dotted 
curve) is able to reproduce the data. Surprisingly, such a LDA [32] was able to reproduce the 
KEK-E307 single-proton data of Ref. [4] with values of T n /r p smaller than 1, as reported in 
the Introduction. This indicates that the two sets of data are inconsistent with each other. 
A clarification of this discrepancy would be desirable. 

In contrast, the single-neutron spectrum measured by KEK-E508 [6] is compatible, 
within error bars, with the previous distribution measured by the KEK-E369 experiment 
[5], as can be seen in Fig. 5. Here, results of the present LDA calculation are also shown, 
both for the primary neutrons coming from the lAMnduced mechanism (dashed line) and 
after including 2AMnduced decays and FSI effects (solid line). The full result is in good 
agreement with both sets of data. 

We note in passing that a model having a sufficiently large value of r n /r p to reproduce 
the KEK-E508 proton spectrum of Fig. 4 would inevitably overestimate the neutron spectra 
of KEK-E508 and KEK-E369 shown in Fig. 5. We would also like to point out that the 
KEK-E508 proton spectrum of Fig. 4 have suffered a much stronger correction from energy 
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FIG. 4. Single-proton kinetic energy spectra for the non-mesonic weak decay of \ 2 C after 
2iV-induced decays and FSI effects are included. Data are from KEK-E508 [6]. All results are 
normalized per non-mesonic weak decay. See detailed explanation in the text. 



— ■ IN without FSI 
— 1N+2N+FSI 
• KEK-E508 
o KEK-E369 




FIG. 5. Single-neutron kinetic energy spectra for the non-mesonic weak decay of ^ C. The 
dashed line corresponds to the distribution of primary neutrons, while the continuous line is ob- 
tained once 2AMnduced decays and FSI effects are included. Data are from KEK-E369 [5] and 
KEK-E508 [6]. The spectrum of primary neutrons (experimental data and the full theoretical 
result) is (are) normalized per liV-induced (total) non-mesonic weak decay. 
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losses in target and detector materials than the neutron spectra of Fig. 5. 

We would like to comment now on a disagreement between theory and experiment con- 
cerning the decay of \B.e. The calculation of Ref. [15] predicted a quite pronounced peak at 
T n ~ 75 MeV in the single-neutron energy spectra from \He non-mesonic decay. Even if a 
LDA formalism is not the best description for this light hypernucleus, we have checked that 
a very similar peak is also produced by the present model. In contrast, almost no peaking 
structure was found in the KEK-E462 experiment of Ref. [6] . 

To gain some insight into this discrepancy, we put forward a simple theoretical argument 
that shows that single-neutron spectra are expected to be less influenced by FSI than single- 
proton spectra. In this discussion we safely neglect 2A r -stimulated decays. We start by 
noting that, for any value of T n /T p , the number of primary neutrons is larger than that of 
primary protons, i.e., N™ d /N p d = 2T n /T p + 1 is always larger than 1. Consequently, due to 
np — * np reactions, the proportion of secondary protons in the proton spectrum N p (mainly 
at low kinetic energies) is larger than the proportion of secondary neutrons in the neutron 
spectrum N n . Note that nn — > nn and pp — > pp reactions, occurring with almost identical 
cross sections, produce the same proportion of secondary neutrons in N n and secondary 
protons in N p , respectively. For this reason, our maxima at T N ~ 75 MeV for ^He are more 
evident for neutrons than for protons (compare Fig. 3 and 4 of Ref. [15]). Due to stronger 
FSI, such maxima completely disappear for \ 2 C (compare Fig. 5 and 6 of Ref. [15]), but 
again protons are more affected by FSI than neutrons. 

According to this discussion, a similar shape of neutron and proton spectra — as indicated 
by the experiment of Ref. [6] for \ 2 C — is only possible for N™ d ~ N p d and thus for small 
values of T n /T p . For T n /T p = we predict N n ~ N p (see Fig. 7), but KEK obtained 
N n /Np ~ 2 for both ^He and \ 2 C. The similar neutron and proton KEK-E462 spectral 
shapes thus look surprising. 

Another point raised in the Conclusion of Ref. [6] has to be commented. There, Okada 
et al. claim that the observed single-neutron spectrum from ^He "indicates the importance 
of the multi-nucleon induced process in the non-mesonic weak decay or/and a large FSI 
effect". Here we want to emphasize that only a very large and quite certainly unrealistic 
proportion of 2A^-stimulated decays could eliminate the maximum that we found at T n ~ 75 
MeV for |He decay. Besides, the disappearance of this maximum would require a too strong 
amount of FSI. 

In Fig. 6 we show the single-neutron energy spectrum obtained for the non-mesonic 
decay of ^ 9 Y. This result is compared with the distribution obtained by the experiment 
KEK-E369 [5]. The quite good agreement between theory and data for a hypernucleus 
as heavy as ^ 9 Y is an indicator of the reliability in using the intranuclear cascade code of 
Ref. [32] to simulate the nucleon FSI in hypernuclear decay. 
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FIG. 6. Single-neutron kinetic energy spectra for the non-mesonic weak decay of ^ 9 Y. Data 
are from KEK-E369 [5]. All results are normalized per non-mesonic weak decay. 



B. Ratio N n /N p 

We would like to compare now our results for the ratio between the number of neu- 
trons and the number of protons produced in the decay of \ 2 C with experimental observa- 
tions. Since our definitive aim is to determine T n /T p by such a comparison, it is convenient 
[15] to start by introducing the number of nucleons of type N (N = n or p) produced in 
n-induced (iV^ Bn ), p-induced (iV^ Bp ), nn-induced (iV^ Bnn ), np-induced (iV^ Bnp ) and pp- 
induced (N^ pp ) decays. By normalizing these quantities per n-, p-, nn-, np- and pp- 
induced decay, respectively, the total number of nucleons of the type N normalized per 
non-mesonic weak decay is given by: 

ATlBn-p I julBpp _|_ Ar2Bnnp , 7\r2Bnp F . ,ir2Bpp F 
>- n \ >- p \ >- nn \ >- np \ pp 

By definition, the nucleon numbers iV^ Bn , iV^ Bp , etc, are independent of the model employed 
to describe the weak decay. Their values depend on the strong interaction part of the 
problem, related to nucleon FSI, and on the framework (finite nucleus or nuclear matter) 
used for treating the hypernuclear structure effects, which produce different phase space 
factors. The dependence on the weak decay model enters Eq. (3.1) via the various partial 
decay widths. 

In table II we report our results for the weak decay model independent nucleon numbers 
in the case of a kinetic energy threshold for nucleon detection of T^ 1 = 60 MeV. From 
Eq. (3.1) and our results of Tables I and II we then determine: 
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TABLE II. Predictions for the weak interaction model independent quantities iVjy , N N p , 
A^ Bnn , iV^ Bnp , A^ Bpp and for Nff of Eqs. (3.1) and (3.4)-(3.6) for fC and nucleon kinetic energies 
T N > 60 MeV. 



ATlBn 


/V 1B P 


jy2Bnn 


/V 2Bn P 




N 2B 


0.91 


0.47 


0.67 


0.50 


0.17 


0.45 


N lBn 


iV p 1Bp 


Ar2Bnn 
1 \p 


jy^Bnp 


N 2Bpp 


N 2B 


0.11 


0.54 


0.08 


0.21 


0.56 


0.27 



^ = 1.33. (3.2) 

We have to note that this result is close to the ones obtained in the finite nucleus calculation 
of Ref. [15], namely N n /N p = 1.38 and 1.42, using the OMEa and OMEf model, respectively. 
If the 2iV-stimulated decay mode is neglected, the present calculation predicts 

(-) =1.28. (3.3) 

These results underestimate the value N n /N p = 2.00 ± 0.17 obtained by KEK-E508 [6] 
for T$ = 60 MeV. Such an occurrence is related to the disagreement between theory and 
experiment for the single-proton spectra of Fig. 4. Our results of Eqs. (3.2) and (3.3) should 
also be compared with the previous experimental determination, N n /N p = 1.73 ± 0.22, 
obtained from KEK-E369 and KEK-E307 data [5] and for a detection threshold of about 
40 MeV. 

It is interesting to mention that, on the contrary, for |He, KEK-E462 [6] measured 
N n /N p = 2.17 ± 0.22, which fairly agrees with the values obtained in Ref. [15], namely 
N n /N p = 1.78 and 1.98 for the OMEa and OMEf models, respectively. 

In order to determine T n /T p for j^C, we now consider a weak decay model independent 
analysis of the mentioned KEK-E508 data. From Eq. (3.1) written for neutrons and protons 
one obtains: 

AriBn^ + ^ibp + M (N™™F nn + N^T np + ivf pp r pp ) 
7V n = l p +i V til _ (3 4) 

NP ATlBnElL + tflBp + _L ^ + ^ (^B nn Fnn + ^Bnpp^ + ^BpppJ 

The ratio T n /T p is thus obtained in terms of our theoretical values for r nn /ri, T np /Yi and 
Tpp/Ti and from data for N n /N p as: 



jV 1B P + n 2B — - ( N 1Bp + N 2B —) — 

1 - ( ATlBn + N 2B^A Nn _ 1B n _ ^B^ 



(3.5) 



p TJ N p n n Y 1 
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where: 



/\7-2Bnnp _i_ A/-2Bnpp . iy-2Bpp-p 
T 2B _ ly N nn ' N n P ' N 1 PP 



* ~ r +r +r ' 1 j 

Using our predictions of Tables I and II together with the datum of KEK-E508 [6], N n /N p = 
2.00 ±0.17, we obtain: 

^ = 0.95 ±0.21. (3.7) 
1 p 

Neglecting the 2iV-stimulated channel the result is: 

= 0.88 ±0.16, (3.8) 
while enhancing arbitrarily the 27V-induced rates by a factor of two we obtain: 

= 1.02 ±0.27. (3.9) 




- r \r 2 ->2r 2 

J- n \ 



The sensitivity of the ratio T n /T p to the values of the 2AMnduced decay widths turns 
out to be moderate, especially if the error bars originated by the datum adopted for N n /N p 
are taken into account. This is due not only to the negligible role of the 2iV-stimulated 
processes in Eq. (3.5) (N^^ and are larger than iV^^/Ti since a quite high energy 

threshold is employed), but also to the particular value of N n /N p used in the analysis, which 
causes a certain cancellation among the 2iV-stimulated contributions in both the numerator 
and denominator of Eq. (3.5): N% B ~ N p B N n /N p . Such a cancellation would be complete, 
thus leading to the extraction of the same central value of T n /T p for any value of T 2 , if a 
ratio N n /N p = N% B /N p B = 1.65 were used in the analysis. 

These occurrences can be understood from Figure 7, which shows the relation between the 
observable ratio N n /N p and T n /T p for different choices of I^/Ti. The dotted line corresponds 
to the case in which 2AMnduced decays and FSI are neglected. Once FSI are incorporated, 
quite different dependencies are obtained. The dot-dashed line refers to the calculation in 
which T 2 is set to 0, the continuous line to the r 2 /Ti ratio predicted by the present LDA 
model and the dashed line to the case in which the size of T 2 is arbitrarily doubled. 

Note in particular that the previously extracted T n /T p values are different from the value 
0.5 given in Ref. [6] and expected on the basis of the relation T n /T p = (N n /N p — l)/2, which 
holds if 2iV-induced decays and FSI effects are ignored, i.e., if A^ Bn = 2, A^ Bp = N p 1Bp = 1, 
N p 1Bn = N^ p f = N^ p = N™ p f = in Eqs. (3.5) and (3.6) (see the point where the 
experimental datum intersects the dotted line in Fig. 7). Contrary to the claim of Ref. [6] 
and according to our results, FSI turn out to be rather important even when one discusses 
the number ratio N n /N p (for which part of the FSI effect is certainly canceled out) and uses 
high kinetic thresholds such as 60 MeV. On the contrary, it is well manifest from Fig. 7 



16 




FIG. 7. Dependence of the observable ratio N n /N p on T n /T p and T2/T1 for a nucleon energy 
threshold of 60 MeV. The horizontal lines show KEK-E508 data [6]. See text for further details. 



that the 2iV-induced decay mechanism plays a relatively small role in the whole range of 
reasonable T n /T p values. 

We also note that a naif estimation of FSI, such as that performed by Kim et al. [5] 
(which also neglect the 2iV-induced mechanism) taking a ratio g/f = 0.11, where / is the 
loss factor of each nucleon type and g the influx factor of one nucleon type converted from 
the other type, gives T n /T p = 0.71±0.14 (see Eq.(4) of Ref. [5]). Although, due to error bars, 
this result is quite compatible with the one obtained here with a more realistic treatment 
of FSI and disregarding 2iV-stimulated decays, r ra /r p = 0.88 ± 0.16, clearly, the difference 
between the two central values demonstrates the following fact: that a detailed account of 
FSI do not support the method used in Ref. [5] to determine Y n /Y p . 

We finalize this subsection by remarking that large ratios such as those of Eqs.(3.7)- 
(3.9) turn out to strongly overestimate the value (0.285) predicted by the weak decay model 
employed in the present paper. Actually, no calculation performed up to date reproduces 
these determinations, which are more in agreement with the previous single-proton data of 
Refs. [4,29,30], as summarized in Table I. Only the presence of an isoscalar, spin-independent 
central operator in the weak transition potential can reproduce larger T n /T p values, as can 
be inferred from the effective treatment of Ref. [26] . 

C. Double coincidence nucleon spectra 

Now we discuss the NN coincidence spectra obtained from our model for the decay of 
a 2 C. Figs. 8 and 9 show, respectively, the distribution of nn and np pairs, as a function 
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FIG. 8. Opening angle distribution of nn pairs normalized per non-mesonic weak decay. 



of the cosine of the opening angle, where FSI effects have been incorporated and a kinetic 
energy cut of 30 MeV has been applied. We observe that the distribution of pairs from 
the liV-induced processes (dash-dotted line) is more back-to-back dominated than that 
coming from the 2AMnduced processes (dashed line). In any case, the liV-induced channel 
still provides the larger contribution of pairs in the whole range of opening angles. 

For cos Onn < —0.4 the results of Figs. 8 and 9 reasonably reproduces the ones of the 
previous finite nucleus calculation reported in Fig. 5 of Ref. [14] and Fig. 9 of Ref. [15]. On 
the contrary, for cos 9 NN > —0.4 a discrepancy is evident, the finite nucleus distributions 
being nearly flat in this region and the nuclear matter ones going monotonously to almost 
vanishing values with cos 9^^. This is ascribable to the different phase space treatment in 
the two approaches. Indeed, the fact that the nuclear matter calculation exhibits a more 
pronounced Fermi motion effect (see Fig. 3) gives rise to a smaller number of outgoing 
nucleons with respect to the finite nucleus case (see Fig. 4). This, in turn, is responsible for 
the fact that the final nucleons of the finite nucleus calculation are more distorted by FSI 
than the nucleons in the nuclear matter approach. As a final result, the NN opening angle 
distributions are more back-to-back correlated in the nuclear matter calculation. 

Our distributions of Figs. 8 and 9 can also be compared with those obtained by KEK- 
E508 and shown in Figure 3 of the recent preprint by Kim et al. [10]. Due to the limited 
statistics of data, we concentrate on the angular region with cos 9nn < —0.7. In this 
region, experiment predicts N nn = 0.083 ± 0.014 and N np = 0.138 ± 0.014, whereas our 
corresponding results are N nn = 0.111 and N np = 0.300, these numbers being normalized 
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FIG. 9. Opening angle distribution of np pairs normalized per non-mesonic weak decay. 



per non-mesonic weak decay. While there is decent agreement in the case of N nn , for N np 
we overestimate the datum by a factor of about 2. The origin of such a discrepancy is 
likely the same as the one responsible for the disagreement on the single-proton spectra of 
Fig. 4. Another indication supporting this hypothesis comes from comparing our result for 
the number of proton-proton pairs for T^ 1 = 30 MeV and cos 0nn < —0.7, N pp = 0.050, 
with the experimental value N pp = 0.005 ± 0.002. One cannot exclude the possibility that 
the experiment systematically underestimated the number of protons emitted in j^C decay, 
thus leading to an underestimation of N p , N np and N pp spectra. On the contrary, for the 
observables involving only neutron detection, N n and N nn , theory and experiment are in 
reasonable accordance. 

The nn and np pair distribution as functions of the total kinetic energy of the pair is 
shown in Figs. 10 and 11, respectively, where the energy of each nucleon is larger than 
a threshold kinetic energy of 30 MeV. The upper panels show the distributions obtained 
without any cut in the opening angle, while in the bottom panels only the back-to-back 
events are kept by applying the restriction cos 9nn < —0.7. We observe that the 2N- 
induced events (dashed lines) in the upper panels scarcely contribute at the position of the 
the primary peak of the lAMnduced contribution (dot-dashed lines). Instead, they enhance 
the total distribution at a pair energy of around 100 MeV and generate a secondary peak 
there, which might become even larger (see the case of the nn pairs) than the primary one 
at the Q-value of about 155 MeV. When the angular cut is applied, many of the events in 
the low energy region are removed and the so-called back-to-back peak at 155 MeV stands 
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FIG. 10. Distribution of the total kinetic energy of nn pairs normalized per non-mesonic weak 
decay. The energy threshold is 30 MeV for each nucleon in the pair. In the upper panel there is 
no angular restriction, while the condition cos 9 nn < —0.7 has been imposed in the results of the 
lower panel. 

out more clearly. However, there is still an important fraction of the events (about 2/3 of nn 
pairs and about 1/2 of np pairs) that lie outside this peak. This is in quantitative agreement 
with the finite nucleus results of our previous works (see Fig. 3 of Ref. [14] and Fig. 10 of 
Ref. [15]), the only qualitative difference being the width of the back-to-back peak, which 
appears more smeared out in the present work due to a more marked effect of Fermi motion. 
In particular, note how the distributions from lV-induced decays of Figs. 10 and 11 extend 
above the Q-value of the non-mesonic decay. 

D. 2N induced strength 

In this subsection we compare the results of our microscopic model for the 2V-induced 
channel with those of the phenomenological model of Ref. [33]. We also analyze here the 
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FIG. 11. Distribution of the total kinetic energy of np pairs normalized per non-mesonic weak 
decay. The energy threshold is 30 MeV for each nucleon in the pair. In the upper panel there is 
no angular restriction, while the condition cos 9 np < —0.7 has been imposed in the results of the 
lower panel. 
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distribution of NN pairs from the various sources of 2iV-induced strength. 

We start by commenting on the claim stated in Ref. [6] about the existence of two 
different types of 2AMnduced decay. One would be a "three-body reaction", in which the 
available energy is equally shared among the three final nucleons, while the other one should 
be due to the absorption of the pion emitted at the weak vertex by a correlated two-nucleon 
pair. According to Ref. [6], in our previous calculations [14,15], where the 2AMnduced 
strength was taken from the model of Ref. [33], we only included the latter. Here we want 
to clarify that there is in fact only one type of 2AMnduced non-mesonic decay and that the 
particular kinematical properties of the emitted particles, prior to FSI, are dictated by the 
dynamics of the 2AMnduced decay mechanism. In the one-pion-exchange model of Ref. [33] 
the absorption by a correlated nucleon pair of the virtual pion emitted at the A vertex is 
driven by a phenomenological two particle-two hole polarization propagator H.2p2h{<f \ q)- In 
the limit (q°, q) — > (m n , 0) this propagator describes the absorption of real pions in nuclear 
matter but we note that it was conveniently extended, via a phase space correction, to 
account for all possible q) values of the exchanged virtual pion. Nevertheless, the pion 
at the A vertex turns out to be preferentially close to its mass shell and, consequently, the 
nucleon emitted at the same vertex has very little kinetic energy left, whereas the other two 
nucleons of the 2AMnduced channel are quite energetic and come out in a back-to-back 
geometry. This is illustrated in Fig. 12, that shows the distribution of momentum values for 
each of the three nucleons emitted by the 2AMnduced mechanism, prior to FSI effects. The 
left panel shows the distribution obtained from the phenomenological model. It is clearly 
seen that the nucleon emitted at the A vertex, whose momentum is denoted by p±, is very 
slow while the other two are equally fast, having momenta on average of about 400 MeV/c. 
In contrast, in the microscopic one-meson-exchange model of Ref. [20] used in the present 
work, the nucleon emitted at the bubble that absorbs the exchanged meson (see Fig. 2) 
can reach high p 2 momentum values by combining the momentum q of the virtual meson 
plus the momentum p' 2 of the correlated nucleon that can be large due to the short range 
nature of the AW interaction. The momentum pi, carried by the nucleon emitted at the 
A vertex, and the momentum p$, corresponding to the nucleon from the spectator bubble, 
will acquire values characteristic of the range of the interaction they come out from, which 
is, respectively, the weak and strong one-meson-exchange models used in Ref. [20]. As 
we observe in the right panel of Fig. 12, the distributions of p\ and p 3 momentum values 
turn out to be quite similar, peaking around 300 MeV/c, while that for p 2 peaks at higher 
momentum values close to 500 MeV/c. 

We now compare the phenomenological and microscopic models after the three primary 
nucleons emitted in the weak decay process are allowed to undergo collisions with the other 
nucleons as they move out of the nucleus. We restrict here to the np-induced decay mode, 
which is the most important one in the microscopic approach and the only one considered by 
the phenomenological model. The distribution of np pairs, normalized per its corresponding 
np-induced transition rate, is shown in Fig. 13 as a function of the opening angle, where 
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the threshold kinetic energy is 30 MeV for each nucleon of the pair. The dotted (solid) line 
shows the results without (with) FSI effects. It is clear that, at variance to the microscopical 
approach, the phenomenological model produces a distinct back-to-back distribution which 
stands quite clearly even when FSI effects are included. The distribution of pairs as a 
function of the total kinetic energy and for a threshold T$ = 30 MeV is displayed in 
Fig. 14. We observe that the amount of N np pairs per np-induced decay event in the absence 
of FSI (area under the dotted curves) is smaller in the phenomenological model, since the 
slow nucleon is always eliminated by the kinetic energy cut of 30 MeV. This situation is 
compensated when the opening angle cut is also applied (dashed lines), since it removes 
more events in the microscopic distribution, which is not so back-to-back dominated. The 
conclusion is that, even if the kinematics of the primary nucleons look quite different, at the 
end, once the effect of FSI is considered and the energy and angular cuts are applied, both 
models produce similar neutron-proton angular and energy spectra per np-induced decay 
event. 

We now discuss the distribution of the various contributions (nn-, np- and pp-induced) 
to the 2iV-induced non-mesonic decay in the microscopic model. The opening angle distri- 
bution of nn and np pairs, again for T$ = 30 MeV, is shown, respectively, on the left and 
right panels of Fig. 15. The angular distribution of both nn and np pairs decreases smoothly 
with increasing cos O^n- We observe that the most important contribution is that of the 
np-induced decay, as we expected from the values of the decay rates shown in Table I. We 
note also that there is a small amount of nn pairs from the pp-induced process and np pairs 
from the nn-induced one which would be zero in the absence of FSI. 

The total pair energy distribution from the 27V-induced channels is shown in Fig. 16, 



23 



phen. 



0.15 



0.10 - 



0.05 



~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



J I I I I I I I I I I I I I I I L 



micr. 

~\ — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

noFSI 
— FSI 



T N > 30 MeV 



j i i i I i i i i I i i i i I LIl 



0.00 

-1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 

cos cos 9 

np np 

FIG. 13. Opening angle distribution of np pairs from the np-induced channel, normalized per 
np-induced decay. 
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FIG. 14. Total kinetic energy distribution of np pairs from the np-induced channel, normalized 
per np-induced decay. 
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FIG. 15. Various 2AMnduced contributions (27V = ran, np,pp) to the opening angle distribution 
of rara (left) and np pairs (right). All results are normalized per total non-mesonic weak decay. 



where an angular cut of cos 9nn < —0.8 has been applied in addition to Tjf = 30 MeV. 
The dominance of the np-induced decay is obvious from figures 15 and 16, although after 
applying the usual kinetic and angular cut-offs the strength of some channels, as e.g. the 
pp-channel contributing to the number of np pairs, can represent up to 20% of the total. 



E. Ratio N nn /N, 



np 



We now discuss a comparison of our predictions for the ratio between the number of nn 
and np pairs for j^C with experimental determinations by KEK-E508 [7-10]. As in previous 
papers [14,15], we introduce the numbers of NN pairs (NN = nn, np or pp) coming from 
one-nucleon induced (Nj^ and N^) and two-nucleon induced (N^ n , N^ p and N^ p ) 
processes, each one of them being normalized per the rate of the corresponding process. 
These weak decay model independent quantities are given in Table III for nucleon kinetic 
energies T N > 30 MeV and two angular regions. The total number of AW pairs emitted per 
non-mesonic weak decay event can be built from: 



N 



NN 



ATlBn-p _|_ /\T 1B Pp _i_ /\T2Bnnp , 7\r2Bnp-p . ,r2Bppp 

ly NN 1 n ' iv NN 1 V ' 1 v NN 1 nn i -< v NN 1 np < NN 1 1 

N]w£n + N^fTp + N N B N T 2 

r n + r p + r 2 ' 



(3.10) 



where: 



ly NN 



A/"2Bnnp 
NN 1 nn 



, yy-2Bnp-p 
' NN L 



np 



7\r2Bpp r 

NN PP 



+ r + r 

~ 1 np ~ 1 pp 



(3.H) 
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FIG. 16. Various 2AMnduced contributions (2N = nn,np,pp) to the pair kinetic energy dis- 
tribution of nn (left) and np pairs (right). All results are normalized per total non-mesonic weak 
decay. 



evidently depends on the weak model employed. 

Before proceeding with the discussion on the ratio N nn /N np , we want to compare the 
results of Table III with the ones obtained within the finite nucleus framework of Refs. 
[14,15] and reported in Table IV. A very good agreement is noticeable for the lAMnduced 
contributions. That is expected and depends on the fact that 1) the considered numbers are 
independent of the weak decay model one uses and 2) the FSI effects are modeled with the 
same intranuclear cascade code in both evaluations. Possible differences have to be ascribed 
to the different frameworks (finite nucleus vs nuclear matter) used for describing hypernu- 
clear structure effects. The disagreement existing for some 2AMnduced contributions, whose 
significance are anyhow relatively low for the total numbers N NN due to the smallness of 
the corresponding decay rates, is due to the different decay channels and (especially) phase 
space factors involved in the two determinations (for a comparison of the phase spaces, we 
refer to the discussion of Fig. 12). 

We now come back to the N nn /N np ratio. From Eq. (3.10) and our results of Tables I 
and III we obtain: 



' nn 

AW 



= 0.36 (3.12) 



for the case with cos 8nn < 0.8. This is in good agreement with the KEK-E508 datum 
0.40 ±0.10 [7]. The comparison of the result (3.12) with those obtained with the calculations 
of Ref. [15] (from Tables I and IV), namely N nn /N np = 0.43 and 0.47, using the OMEa and 
OMEf models, respectively, show the differences in the partial decay rates and phase spaces 
predicted by the different models. If the 2A r -stimulated decay mode is neglected, the present 
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TABLE III. Predictions for the weak interaction model independent quantities Njf£, N N ^, 

N NN n , N NN P > n nn P and for n nn of E q s - ( 3 - 10 ) and ( 3 - u ) for A C > integrated over all opening 
angles and for nucleon energies Tjy > 30 MeV. The numbers in parentheses correspond to the 



angular region with cos 9nn < —0.8. 



ATlBn 
A y nn 


7V 1B P 

±y nn 


;y2Bnn 
Jv nn 


Af2Bnp 

nn 


/\T 2B PP 

nn 


N 2B 
±y nn 


0.55 (0.31) 


0.11 (0.03) 


0.53 (0.15) 


0.24 (0.06) 


0.05 (0.01) 


0.22 (0.06) 


ATlBn 
A 'np 


ATiBp 

1 ' np 


Ar2Bnn 
1 v np 


AT 2Bn P 
A 'np 


/V 2B PP 

1 ' np 


1 ' np 


0.34 (0.09) 


0.64 (0.34) 


0.27 (0.05) 


0.45 (0.10) 


0.40 (0.09) 


0.43 (0.09) 


ATlBn 
pp 


ATiBp 
pp 


]y~2Bnn 
pp 


/\T 2B np 
PP 


/V 2B pp 

pp 


N 2B 

pp 


0.04 (0.01) 


0.15 (0.05) 


0.02 (0.004) 


0.08 (0.02) 


0.30 (0.08) 


0.12 (0.03) 



TABLE IV. Predictions of the finite nucleus calculation of Ref. [15] for the weak interaction 
model independent quantities Nj^, N^j^ and for N$P N (integrated over all angles and for nucleon 
energies T/v > 30 MeV) for ]^C. The numbers in parentheses correspond to the angular region with 
cos6 NN < -0.8. 



ATlBn 

nn 


AT 1B P 

nn 


N 2B 


0.57 (0.31) 


0.11 (0.03) 


0.30 (0.12) 


ATlBn 
np 


ATiBp 

np 


N 2B 

np 


0.34 (0.09) 


0.68 (0.32) 


0.39 (0.10) 


ATlBn 

pp 


7V 1B P 
pp 


N 2B 

pp 


0.04 (0.01) 


0.17 (0.05) 


0.05 (0.01) 
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calculation predicts: 

(N \ 1N 

\jt P ) = °- 34 ' (3,13) 

thus emphasizing a relatively small effect of the 2iV-induced channels on correlation ob- 
servables appropriately chosen. The above result (3.13) should be compared with the naif 
estimation N nn /N np = T n /T p = 0.285 obtained by neglecting both FSI and the 2AMnduced 
decay channel: contrary to what occurs to single-nucleon spectra, i.e., to N n /N p , the ef- 
fect of FSI in the present case of coincidence observables turns out to be not too marked. 
The reliability of the above naif equality improves by using higher thresholds T$ and 
more restrictive back-to-back constraints. We can therefore conclude that, for the cases 
of experimental interest such as those considered here and in Section III B, the equation 
N nn /N np = T n /T p turns out to have a wider range of approximate validity than the assump- 
tion N n /N p = (2r n + T p )/T p sometimes used in experimental analyses, whose utilization to 
estimate r n /r p must be avoided. 

We now consider a weak decay model independent analysis of the N nn /N np data of 
Table V to determine r n /r p . These data have been obtained by KEK for |He [9] and \ 2 C 
[7,10] with a threshold T$ = 30 MeV and various angular restrictions. The ratio r„/r p is 
determined by using theoretical values for T nn /Ti, T np /Ti and r pp /Ti and data for N nn /N np 
from: 

-^2 ( ^lBp , 7\r2B-T2^ N nn 

N 

np (3.14) 

J- 2 

np ' ^'rap-p / Af - 1 ' ran * ' nn -p 

1 1 / ^np J- 1 

where Nf^ N is given by Eq. (3.11). From our predictions of Tables I and III and the ^ 2 C 
datum of KEK-E508 derived with the restrictions T N > 30 MeV and cos 9 NN < —0.8 we 
obtain: 

^ = 0.34 ±0.15. (3.15) 
1 P 

Besides, neglecting the 2iV-stimulated channel: 

l-^j =0.37 ±0.14, (3.16) 

while enhancing arbitrarily the 2iV-induced rates by a factor of two: 

r 2 ->2r 2 

= 0.32 ±0.16. (3.17) 

The ratio (3.15) [(3.16)] must be compared with T n /r p = 0.27 ± 0.14 [r n /r p = 0.36 ± 0.14] 
obtained from the results of Ref. [15] (listed in Table IV) together with a value of r 2 /Ti = 



iV 1B P + N 2B — - ( N 1Bp + N 2B — 

Y nn 1 ran -p yj np 1 np p 

T p / ^IBn , /v2b£2\ ^nn _ ^iBn _ at2b£2 
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TABLE V. iHe [9] and i 2 C [7,10] KEK data for the ratio N nn /N np corresponding to a nucleon 
threshold T$ = 30 MeV and different opening angle regions. The corresponding determinations 
of T n /r p obtained in the present work and from the results of Ref. [15] are also reported. 





IHe 






A u 


Angular range 




r n /T p [15] 




r n /r p [this work] 


COS0NN < —0.9 






0.45 ±0.12 


0.43 ±0.17 


COS0NN < —0.8 


0.45 ±0.11 


0.27 ±0.11 


0.40 ±0.10 


0.34 ±0.15 


COS0NN < —0.7 






0.60 ±0.12 


0.66 ± 0.24 



0.25 used there. While the ratios are very similar when T 2 = 0, a certain difference is 
obtained for those including 2A^-stimulated decay effects. This is quite obvious and signals 
the different models used in the two works for the 2AMnduced channels. 

In Figure 17 we show the relation between N nn /N np and T n /T p for different choices of 
r 2 /Ti and for T$ = 30 MeV and cos6nn < —0.8. The dotted line corresponds to the case 
in which 2AMnduced decays and FSI are neglected. Different behaviours are obtained once 
FSI are incorporated. The dot-dashed line refers to the calculation in which T 2 is set to 0, 
the continuous line to the T2/T1 ratio predicted by the present LDA model and the dashed 
line to the case in which the size of T 2 is arbitrarily doubled. The comparison of Fig. 17 
with Fig. 7 clearly illustrates the fact that FSI affect much more the extraction of r n /r p 
from N n /N p than from N nn /N np . 

Table V summarizes the values of T n /T p derived in our analysis of KEK data for N nn /N np 
in i 2 C. Also reported is the result for iHe of the finite nucleus calculation of Ref. [15]. In 
all cases the values reported for r„/r p have been obtained by including the 2A^-stimulated 
contributions as predicted by the two different approaches. As a signal of the limited statis- 
tics of the data employed, the central value for r n /r p in \ 2 C shows a certain undesired 
dependence on the opening angle region. Nevertheless, the error bars are so big that values 
in the interval 0.4-0.5 are common to all determinations. The final result can thus be given 
by the weighted average of the three partial results: 

best 

= 0.43 ±0.10. (3.18) 




When 2AMnduced decay channels are neglected, a similar analysis of data supplies a very 
close "best" value: 

IN best 

= 0.46 ±0.09. (3.19) 




We want to note that for the case with cosO^n < —0.7, Kim et al. [10] have recently 
extracted a ratio r„/T p = 0.51±0.13±0.04by employing the number of detected proton- 
proton pairs in addition to the measurement of N nn and N np . A schematic method for 
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FIG. 17. Dependence of the observable ratio N n /N p on T n /T p and ^/Ti for a nucleon energy 
threshold of 30 MeV and cos NN < -0.8. The horizontal lines show KEK-E508 data [7]. See text 
for further details. 

accounting for FSI effects, similar to the one employed in Ref. [5] for single-nucleon ob- 
servables, was adopted for this purpose. In addition, 2iV-stimulated decay channels were 
neglected. The discrepancy between this KEK determination and our detailed prediction of 
Table V for cos# < —0.7 (compare especially the central values) demonstrate that analyses 
of data such as those of Refs. [5,10] show certain limitations. 



We have presented a study of single- and double-coincidence nucleon spectra for the 
non-mesonic weak decay of A-hypernuclei. One-meson-exchange models have been used 
to describe one- and two-nucleon induced decay processes, AN — > nN and ANN — > nNN, 
in a nuclear matter framework which has been adapted for finite nuclei predictions via 
the local density approximation. Besides, an intranuclear cascade code based on Monte 
Carlo techniques has simulated the final state interactions of the outgoing nucleons with the 
residual nucleus. 

Unlike previous papers, here we have adopted a microscopic approach for describing the 
two-nucleon induced decay, which included the channels Ann — > nnn and App — > npp in 
addition to the mode Anp — > nnp already considered in previous phenomenological studies. 
The use of a different theoretical approach for calculation meets our main purpose of making 
the extraction of the ratio r n /r p from data less model dependent. 

The results have been compared with previous finite nucleus analyses [14,15], with special 
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consideration for the two-nucleon induced decay mechanism, and with a considerable amount 
of recent data by KEK [5-10]. Apart from some difference ascribable to the phase space 
dependence and to the different weak decay models adopted for the two-nucleon induced 
decay channels, the present predictions for the observable ratios N n /N p and N nn /N np confirm 
the finite nucleus results. 

The single-neutron spectra for \ 2 C and fY measured by KEK-E369 and KEK-E508 have 
been reproduced with reasonable accuracy. On the contrary, the KEK-E508 single-proton 
spectrum for ^ 2 C is in strong disagreement with our prediction. The theoretical distribution 
would reproduce these data only by enforcing artificially large values for T n /T p , which, in 
turn, would lead to a serious overestimation of the mentioned single-neutron spectra. As 
a consequence of that, through the weak interaction model independent analysis of N n /N p 
data for \ 2 C we determined T n /T p = 0.95 ± 0.21, a value that largely overestimates pure 
theoretical predictions ranging between 0.3 and 0.5. 

The contribution of the two-nucleon induced decay channels in these analyses turn out 
to be of moderate size when high nucleon kinetic energy cuts such as 60 MeV are imposed. 
On the contrary, FSI effects reveal to be of great importance when determining T n /T p 
from single-nucleon data. A demonstration of that is the inapplicability of the relation 
T n /T p = (N n /N p — l)/2, which neglects FSI and two-nucleon stimulated decay effects, 
sometimes used in experimental analyses. Even the naif model of FSI proposed by Kim et 
al. in Ref. [5] is far from being supported by a detailed calculation. 

Concerning nucleon-nucleon correlation observables, while there is fair agreement among 
our predictions and data in the case of N nn , we overestimate the observations in the case 
of N np . The origin of such a difference is likely the same as the one responsible for the 
mentioned disagreement on the single-proton spectra. Another indication supporting this 
possibility comes from the fact that we also overestimate the experimental distributions for 
N pp . We have thus advanced the hypothesis that the experiment could have systematically 
underestimated the number of protons emitted in \ 2 C decay, thus leading to the underes- 
timation of N p , N np and Npp spectra. On the contrary, for the observables involving only 
neutron detection, N n and N nn , theory and experiment turn out to agree. The same con- 
clusion can be drawn in terms of the finite nucleus results. A clarification of the mentioned 
eventuality is desirable for an accurate determination of the T n /T p ratio. 

Indeed, in this paper as well as in previous ones [14,15], r„/r p has been determined from 
data on N nn /N np . Here, for \ 2 C we have obtained T n /T p = 0.43 ± 0.10, whose central value 
could be lowered to about 0.3 if the measured neutron-proton distributions resembled those 
of our evaluations. The conclusion remains the same by using the finite nucleus analysis. 
The effect of the two-nucleon induced decay modes on the extracted T n /T p turn out to 
be small due to the restrictions imposed on angular and energy correlations. In any case, 
one has to bear in mind that ratios in the interval 0.3-0.5 would be compatible with most 
of today's pure theoretical evaluations, which are also affected by uncertainties, especially 
due to the degree of arbitrariness in the experimentally unknown baryon couplings. These 
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models can also reproduce the observed total non-mesonic decay rates. 

Before concluding, we want to make a brief comment on the two-nucleon induced mode. 
Because of the small or anyhow moderate dependence of the determined T n /T p values on 
these decay channels, triple coincidence measurements reveal to be necessary for the purpose 
of disentangling one- and two-nucleon stimulated decays in observed events. 

To summarize, while we can safely assert that analyses of correlation measurements 
definitely solve the longstanding puzzle on the ratio T n /T p , single-nucleon spectra studies 
still provide ratios which are incompatible with what is obtained from pure theoretical 
models. We hope it will be possible to clarify soon the reasons of such a discrepancy. 

ACKNOWLEDGMENTS 

This work is partly supported by EURIDICE HPRN-CT-2002-00311, contract FIS2005- 
03142 from MEC (Spain) and FEDER, INFN-MEC collaboration agreement number 06-36, 
and Generalitat de Catalunya contract 2005SGR-00343. This research is part of the EU 
Integrated Infrastructure Initiative Hadron Physics Project under contract number RII3- 
CT-2004-506078. 



32 



REFERENCES 



[1] W. M. Alberico and G. Garbarino, Phys. Rep. 369, 1 (2002); in Hadron Physics, IOS 
Press, Amsterdam, 2005, p. 125. Edited by T. Bressani, A. Filippi and U. Wiedner. 
Proceedings of the International School of Physics "Enrico Fermi", Course CLVIII, 
Varenna (Italy), June 22 - July 2, 2004. 

[2] E. Oset and A. Ramos, Prog. Part. Nucl. Phys. 41, 191-253 (1998). 

[3] O. Hashimoto et al, Phys. Rev. Lett. 88, 042503 (2002). 

[4] Y. Sato et al., Phys. Rev. C 71, 025203 (2005). 

[5] J. H. Kim el al., Phys. Rev. C 68, 065201 (2003). 

[6] S. Okada et al, Phys. Lett. B 597, 249 (2004). 

[7] H. Bhang, in DAPHNE 2004: Physics at meson factories, Frascati Phys. Ser. 36, 243 
(2005). Edited by F. Anulli, M. Bertani, G. Capon, C. Curceanu-Petrascu, F. L. Fabbri 
and S. Miscetti. 

[8] H. Outa, in Hadron Physics (Ref. [1]) p. 219; H. Outa et al, Nucl. Phys. A 754, 157c 
(2005). 

[9] B. H. Kang et al., Phys. Rev. Lett. 96, 062301 (2006). 
[10] M. J. Kim et al., nucl-ex/0601029. 

[11] M. Agnello et al., Nucl. Phys. A 754, 399c (2005); T. Bressani, E. Botta, A. Feliciello 

and V. Paticchio, Nucl. Phys. A 754, 410c (2005). 
[12] T. Nagae, Nucl. Phys. A 754, 443c (2005). 

[13] I. Tanihata et al., Hypernuclei with Stable Heavy Ion Beam and Rl-beam Induced 

Reactions at GSI (HypHI), Letter of Intent, February 2005. 
[14] G. Garbarino, A. Parreho and A. Ramos, Phys. Rev. Lett. 91, 112501 (2003). 
[15] G. Garbarino, A. Parreho and A. Ramos, Phys. Rev. C 69, 054603 (2004); W. M. 

Alberico, G. Garbarino, A. Parreho and A. Ramos, nucl-th/0407046, in DAPHNE 2004: 

Physics at meson factories, Frascati Phys. Ser. 36, 249 (2005). Edited by F. Anulli, M. 

Bertani, G. Capon, C. Curceanu-Petrascu, F. L. Fabbri and S. Miscetti. 
[16] D. Jido, E. Oset and J. E. Palomar, Nucl. Phys. A 694, 525 (2001). 
[17] A. Parreho and A. Ramos, Phys. Rev. C 65, 015204 (2002); A. Parreho, A. Ramos and 

C. Bennhold, Phys. Rev. C 56, 339 (1997). 
[18] K. Itonaga, T. Ueda and T. Motoba, Phys. Rev. C 65, 034617 (2002). 
[19] E. Bauer and F. Krmpotic, Nucl. Phys. A 717, 217 (2003). 
[20] E. Bauer and F. Krmpotic, Nucl. Phys. A 739, 109 (2004). 
[21] S. Ajimura et al., Phys. Rev. Lett. 84, 4052 (2000). 
[22] T. Maruta et al, nucl-ex/0509016; Nucl. Phys. A 754, 168c (2005). 
[23] W.M. Alberico, G. Garbarino, A. Parreho and A. Ramos, Phys. Rev. Lett. 94, 082501 

(2005). 

[24] C. Barbero, A. P. Galeao and F. Krmpotic, Phys. Rev. C 72, 035210 (2005); C. Barbero, 
C. Dc Conti, A. P. Galeao and F. Krmpotic, Nucl. Phys. A 726, 267 (2003). 



33 



[25] K. Sasaki, M. Izaki, and M. Oka, Phys. Rev. C 71, 035502 (2005); K. Sasaki, T. Inoue 

and M. Oka, Nucl. Phys. A 707, 477 (2002). 
[26] A. Parreno, C. Bennhold and B.R. Holstein, Phys. Rev. C 70, 051601(R) (2004); Nucl. 

Phys. A 754, 127c (2005). 
[27] A. Montwill et al, Nucl. Phys. A 234, 413 (1974). 
[28] A. Sakaguchi et al., Phys. Rev. C 43 (1991) 73. 
[29] J. J. Szymanski et al., Phys. Rev. C 43, 849 (1991). 
[30] H. Noumi et al., Phys. Rev. C 52, 2936 (1995). 

[31] H. Noumi et al., Proceedings of the IV International Symposium on Weak and Electro- 
magnetic Interactions in Nuclei, edited by H. Ejiri, T. Kishimoto and T. Sato (World 
Scientific, Singapore, 1995), p. 550. 

[32] A. Ramos, M. J. Vicente- Vacas and E. Oset, Phys. Rev. C 55, 735 (1997); 66, 039903(E) 
(2002). 

[33] A. Ramos, E. Oset and L. L. Salcedo, Phys. Rev. C 50, 2314 (1994). 
[34] W. M. Alberico, A. De Pace, G. Garbarino and A. Ramos, Phys. Rev. C 61, 044314 
(2000). 

[35] E. Oset and L. L. Salcedo, Nucl. Phys. A 443, 704 (1985). 

[36] W. M. Alberico, A. De Pace, G. Garbarino, R. Cenni, Nucl. Phys. A 668, 113 (2000). 
[37] E. Bauer, A. Polls and A. Ramos, Phys. Rev. C 58, 1052 (1998). 
[38] W. M. Alberico, M. Ericson and A. Molinari, Ann. Phys. C 154, 356 (1984). 
[39] R. Machleidt, K. Holinde and Ch. Elster; Phys. Rep. 149, 1 (1987). 
[40] M. B. Barbara, A. De Pace, T. W. Donnelly and A. Molinari, Nucl. Phys. A 596, 553 
(1996). 

[41] J. Cugnon and P. Henrotte, Eur. Phys. J. A 16, 393 (2003). 



34 



